ABSTRACT

Natural ecosystems across the globe are increasingly suffering from the pressures of on-going global change. Pollinators, in particular, have been reported to decline in abundance and diversity in different parts of the world, with potential implications for many plant species given the dependence of 87.5% of plant species on pollinators for reproduction. Although a strong body of research has focused on pollinator visitation frequency and plant reproduction, most of these efforts have been done at the individual species level, largely ignoring community-level implications. Indeed, interaction network structure can affect pollinator functional roles, their visitation patterns to different plant species and ultimately the pollination service they deliver to plants. Thus, a major knowledge gap at the moment is a mechanistic understanding of how pollinator visitation and the structure of plant-pollinator visitation networks affect plant reproduction. Here, we use 16 well-resolved plant-pollinator networks from Mediterranean shrub ecosystems to evaluate how pollinator species diversity, community composition and interaction network structure affect plant reproduction for 19 species of plants. Further, we link pollinator visitation frequencies with the efficiency of individual pollinator species. To this end, we use data on fruit and seed set for all plant species as well as data on the efficiency of a single visit by different pollinator species on 12 of these plant species. This work represents one of the first efforts linking pollination visitation and plant reproduction from a community-wide perspective using a well-replicated dataset.

INTRODUCTION

POllinators provide key services to plants by facilitating pollen flow between individuals. Given recent declining trends found for some pollinator species in some regions of the planet (Potts et al. 2010), efforts have been devoted to study the effect of decreasing pollinator diversity and abundance for plant reproduction. Many of these approaches focus on pairwise plant-pollinator interactions [doi:10.1111/1365-2745.13055]. However, community-level analyses remain to be conducted. In particular, we have some knowledge on the effect of pollinator species richness for plant reproduction for one species of plant (Albrecht et al. 2012). We also have an understanding of the effects of pollinator diversity and niche complementarity between pollinators for plant reproductive success within an artificially-created community of 9 plant species (Fründ et al. 2013). Frund et al (2013) found that greater species diversity and functional complementarity led to greater average reproduction for the plant species surveyed. However, these approaches are based on experimental research, where plant and pollinator communities are artificially recreated within enclosures and thus may not reflect real-world conditions. In addition, other aspects of community structure and trade-offs between them remain to be explored. Indeed, we know that the diversity of pollinators visiting a plant (the degree in network terms) and the total number of visits they provide, have an important effect for reproductive success. However, we have yet to know whether these variables contribute more to plant reproductive success than others related to the position of a particular plant species within the wider plant community that hosts it. In particular, our study focuses on the importance of pollinator diversity and total number of visits for plant reproductive success in comparison to that of the plant´s position in the community. Specifically, niche overlap between plant species (which could give us a proxy for the existence of competition for pollination between plant species) and the position of the plant given by the connections it establishes with other species, e.g. its centrality within the network, where central species are those receiving the greatest number of connections from other species, could have an influence in a plant´s reproductive success. Furthermore, we do not know how the overall structure of a community, and the architecture of the network of interactions it hosts, affects the reproductive success of the whole community. In particular, one of the networks metrics that has received the most attention in the community ecology literature is that of nestedness. Nestedness is the property by which specialists interact with a subset of the species that generalists interact with. Although there is ongoing debate in the literature, some studies have found that more nested networks are more stable and resilient to perturbations because nestedness implies a certain level of redundancy and thus the possibility for functional replacement in the case of species loss. However, plant reproductive success depends on the delivery of conspecific pollen and thus of a certain level of specialization or niche divergence. Yet at present we do not know how the interplay between redundancy and complementarity determines reproductive success. Finally, average values of reproductive success at the community level can be driven to a large degree by one single plant species. Yet, what will determine the persistence of a diverse plant community is the presence of some sort of “equity” in reproductive success across the whole community. We thus also focus on the effect of community structure for equity in reproductive success across the plant community. Our study aims to answer the following questions: 1) How do “pollinator visits” and “plant role” within the community affect its reproductive success? COMO LLAMAMOS A ESTO??

METHODS

Our study was conducted in SE Spain within the area of influence of Doñana National Park?? Here, we surveyed 16 Mediterranean woodland patches. Each site was surveyed 7 times following a 100-m x 2 m transect for 30 mins. Along each transect we identified all plant species and the pollinators that legitimately visited them. In addition, within each site we evaluated fruit set, seed number per fruit and seed weight for several individuals belonging to 15 plant species. Not all plant species were present in all transects.

DATA ANALYSES

First, we checked for sampling completeness for the pollinator and the plant community as well as for plant-pollinator links using Chao 1 asymptotic species richness estimators. We first estimated the richness of pollinators, plants and plant–pollinator links accumulated as sampling effort increased up to 100% sampling coverage using package iNEXT.

Species-level

The data for all rounds was pooled to construct one plant-pollinator interaction network per site, for which we calculated a set of relevant network metrics. At the species level, we focused on: degree (i.e., number of pollinator species visiting a particular plant), total number of visits, centrality and average niche overlap between a focal species and other plant species in the community. We fit three models, one including only the diversity of pollinator species and the total number of species. One including only metrics related to a plant´s position in the community and a third one combining both sets of models and used information criterion methods to find the best model. We included plant species nested within site as random effect to account for non-independence of several individuals measured for the same plant species.

Network level

We scaled this at the community following the same logic and evaluated the importance of pollinator diversity and visits compared to that of network structure. Including in this case nestedness and functional complementarity as metrics. In the case of nestedness. Many approaches to measure it use null-model approaches to compare across different networks. However, this approach has the problem that z-scores depend strongly on network size and connectance as shown by this recent paper by Song et al. What these authors propose is to use a normalized value of nestedness, and to control by the effect of connectance and network size, thus obtaining a new value of nestedness NODFc which is independent of network size and thus comparable across different networks.

Equity in reproductive success

Finally, we used the same framework to evaluate the contribution of the different variables we consider to equity in plant reproductive success across the community, by including as a response variable the number of plant species that produce 50% of fruit set recorded for our study area at each site.

RESULTS

Within our transects we recorded 57 species of plants and 277 species of pollinators (HAY QUE MIRAR MORFOESPECIES!). Our sampling completeness analyses revealed that with our survey we were able to capture 18-62% of pollinator species (average= 35%), 47-98% for plant species (average= 78%) and 13-41% for plant-pollinator links (average=27%), in line with that obtained with other studies (e.g., Chacoff et al), yet slghtly smaller in the case of pollinators given the great diversity of pollinators found in the Mediterranean region in particular and in the our study area, a hotspot of insect diversity (?¿?).

Species-level

## NULL

##  Family: gaussian  ( identity )
## Formula:          
## tot.fruitset ~ norm_degree + tot.visits.sc + (1 | Site_ID) +  
##     (1 | Site_ID)
## Data: f2
## 
##      AIC      BIC   logLik deviance df.resid 
##    185.0    208.5    -86.5    173.0      366 
## 
## Random effects:
## 
## Conditional model:
##  Groups    Name        Variance Std.Dev.
##  Site_ID   (Intercept) 0.002385 0.04884 
##  Site_ID.1 (Intercept) 0.002385 0.04884 
##  Residual              0.090209 0.30035 
## Number of obs: 372, groups:  Site_ID, 16
## 
## Dispersion estimate for gaussian family (sigma^2): 0.0902 
## 
## Conditional model:
##               Estimate Std. Error z value Pr(>|z|)    
## (Intercept)    0.68325    0.03729  18.322   <2e-16 ***
## norm_degree    0.36912    0.15687   2.353   0.0186 *  
## tot.visits.sc  0.06409    0.10745   0.597   0.5508    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##  Family: gaussian  ( identity )
## Formula:          
## tot.fruitset ~ weigh_closeness + niche.overlap + (1 | Site_ID)
## Data: f2
## 
##      AIC      BIC   logLik deviance df.resid 
##    182.6    202.2    -86.3    172.6      364 
## 
## Random effects:
## 
## Conditional model:
##  Groups   Name        Variance Std.Dev.
##  Site_ID  (Intercept) 0.005978 0.07731 
##  Residual             0.089986 0.29998 
## Number of obs: 369, groups:  Site_ID, 16
## 
## Dispersion estimate for gaussian family (sigma^2): 0.09 
## 
## Conditional model:
##                 Estimate Std. Error z value Pr(>|z|)    
## (Intercept)       0.6482     0.0539  12.026   <2e-16 ***
## weigh_closeness   1.2291     0.5402   2.275   0.0229 *  
## niche.overlap     0.0741     0.2190   0.338   0.7351    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##  Family: gaussian  ( identity )
## Formula:          
## tot.fruitset ~ norm_degree + tot.visits.sc + weigh_closeness +  
##     niche.overlap + (1 | Site_ID)
## Data: f2
## 
##      AIC      BIC   logLik deviance df.resid 
##    183.2    210.6    -84.6    169.2      362 
## 
## Random effects:
## 
## Conditional model:
##  Groups   Name        Variance Std.Dev.
##  Site_ID  (Intercept) 0.004664 0.0683  
##  Residual             0.089661 0.2994  
## Number of obs: 369, groups:  Site_ID, 16
## 
## Dispersion estimate for gaussian family (sigma^2): 0.0897 
## 
## Conditional model:
##                  Estimate Std. Error z value Pr(>|z|)    
## (Intercept)      0.606393   0.057499  10.546   <2e-16 ***
## norm_degree      0.322639   0.176341   1.830   0.0673 .  
## tot.visits.sc   -0.009229   0.118653  -0.078   0.9380    
## weigh_closeness  0.763344   0.566617   1.347   0.1779    
## niche.overlap    0.193543   0.231052   0.838   0.4022    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##    dAIC df
## m2 0.0  5 
## m3 0.6  7 
## m1 2.3  6

##  Family: poisson  ( log )
## Formula:          
## mean.seedn.sc ~ norm_degree + tot.visits.sc + (1 | Site_ID) +  
##     (1 | Site_ID)
## Data: s2
## 
##      AIC      BIC   logLik deviance df.resid 
##    139.8    158.0    -64.9    129.8      276 
## 
## Random effects:
## 
## Conditional model:
##  Groups    Name        Variance  Std.Dev. 
##  Site_ID   (Intercept) 1.544e-09 3.929e-05
##  Site_ID.1 (Intercept) 1.544e-09 3.929e-05
## Number of obs: 281, groups:  Site_ID, 16
## 
## Conditional model:
##               Estimate Std. Error z value Pr(>|z|)    
## (Intercept)    -2.8580     0.4478  -6.383 1.74e-10 ***
## norm_degree     0.2258     2.0847   0.108    0.914    
## tot.visits.sc   0.8924     1.1378   0.784    0.433    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##  Family: poisson  ( log )
## Formula:          
## mean.seedn.sc ~ weigh_closeness + niche.overlap + (1 | Site_ID)
## Data: s2
## 
##      AIC      BIC   logLik deviance df.resid 
##    136.2    150.7    -64.1    128.2      274 
## 
## Random effects:
## 
## Conditional model:
##  Groups  Name        Variance  Std.Dev. 
##  Site_ID (Intercept) 1.643e-09 4.053e-05
## Number of obs: 278, groups:  Site_ID, 16
## 
## Conditional model:
##                 Estimate Std. Error z value Pr(>|z|)    
## (Intercept)      -2.9672     0.5978  -4.963 6.93e-07 ***
## weigh_closeness  -3.1989     6.1777  -0.518    0.605    
## niche.overlap     2.7723     2.0688   1.340    0.180    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##  Family: poisson  ( log )
## Formula:          
## mean.seedn.sc ~ norm_degree + tot.visits.sc + weigh_closeness +  
##     niche.overlap + (1 | Site_ID)
## Data: s2
## 
##      AIC      BIC   logLik deviance df.resid 
##    139.9    161.6    -63.9    127.9      272 
## 
## Random effects:
## 
## Conditional model:
##  Groups  Name        Variance  Std.Dev. 
##  Site_ID (Intercept) 1.631e-09 4.038e-05
## Number of obs: 278, groups:  Site_ID, 16
## 
## Conditional model:
##                 Estimate Std. Error z value Pr(>|z|)    
## (Intercept)      -3.1597     0.7188  -4.396  1.1e-05 ***
## norm_degree       1.2924     2.2634   0.571    0.568    
## tot.visits.sc    -0.0259     1.4097  -0.018    0.985    
## weigh_closeness  -4.5053     6.7168  -0.671    0.502    
## niche.overlap     3.0854     2.5918   1.190    0.234    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##         dAIC df
## m2.seed 0.0  4 
## m1.seed 3.6  5 
## m3.seed 3.6  6

##  Family: poisson  ( log )
## Formula:          
## mean.seedn.sc ~ norm_degree + tot.visits.sc + (1 | Site_ID) +  
##     (1 | Site_ID)
## Data: s3
## 
##      AIC      BIC   logLik deviance df.resid 
##     49.3     67.0    -19.7     39.3      246 
## 
## Random effects:
## 
## Conditional model:
##  Groups    Name        Variance  Std.Dev. 
##  Site_ID   (Intercept) 6.059e-10 2.462e-05
##  Site_ID.1 (Intercept) 6.059e-10 2.462e-05
## Number of obs: 251, groups:  Site_ID, 16
## 
## Conditional model:
##               Estimate Std. Error z value Pr(>|z|)    
## (Intercept)    -4.5458     0.9673  -4.699 2.61e-06 ***
## norm_degree     2.3525     3.7319   0.630    0.528    
## tot.visits.sc   0.3117     3.4548   0.090    0.928    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##  Family: poisson  ( log )
## Formula:          
## mean.seedn.sc ~ weigh_closeness + niche.overlap + (1 | Site_ID)
## Data: s3
## 
##      AIC      BIC   logLik deviance df.resid 
##     47.3     61.4    -19.7     39.3      244 
## 
## Random effects:
## 
## Conditional model:
##  Groups  Name        Variance  Std.Dev. 
##  Site_ID (Intercept) 1.478e-09 3.845e-05
## Number of obs: 248, groups:  Site_ID, 16
## 
## Conditional model:
##                 Estimate Std. Error z value Pr(>|z|)    
## (Intercept)       -3.730      1.085  -3.437 0.000589 ***
## weigh_closeness    1.392     11.579   0.120 0.904274    
## niche.overlap     -2.553      5.021  -0.508 0.611170    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##  Family: poisson  ( log )
## Formula:          
## mean.seedn.sc ~ norm_degree + tot.visits.sc + weigh_closeness +  
##     niche.overlap + (1 | Site_ID)
## Data: s3
## 
##      AIC      BIC   logLik deviance df.resid 
##     50.8     71.9    -19.4     38.8      242 
## 
## Random effects:
## 
## Conditional model:
##  Groups  Name        Variance  Std.Dev. 
##  Site_ID (Intercept) 1.331e-09 3.649e-05
## Number of obs: 248, groups:  Site_ID, 16
## 
## Conditional model:
##                 Estimate Std. Error z value Pr(>|z|)   
## (Intercept)       -3.953      1.304  -3.030  0.00244 **
## norm_degree        1.862      4.014   0.464  0.64271   
## tot.visits.sc      1.956      4.605   0.425  0.67100   
## weigh_closeness   -2.045     13.409  -0.152  0.87881   
## niche.overlap     -3.490      6.164  -0.566  0.57124   
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##          dAIC df
## m2.seedb 0.0  4 
## m1.seedb 2.0  5 
## m3.seedb 3.5  6

##  Family: poisson  ( log )
## Formula:          
## mean.seedweight ~ norm_degree + tot.visits.sc + (1 | Site_ID)
## Data: s2
## 
##      AIC      BIC   logLik deviance df.resid 
##     61.8     76.4    -26.9     53.8      277 
## 
## Random effects:
## 
## Conditional model:
##  Groups  Name        Variance Std.Dev. 
##  Site_ID (Intercept) 1.14e-09 3.376e-05
## Number of obs: 281, groups:  Site_ID, 16
## 
## Conditional model:
##               Estimate Std. Error z value Pr(>|z|)    
## (Intercept)    -4.1003     0.7878  -5.205 1.94e-07 ***
## norm_degree     1.2074     3.4231   0.353    0.724    
## tot.visits.sc   0.5444     2.1169   0.257    0.797    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##  Family: poisson  ( log )
## Formula:          
## mean.seedweight ~ weigh_closeness + niche.overlap + (1 | Site_ID)
## Data: s2
## 
##      AIC      BIC   logLik deviance df.resid 
##     61.7     76.2    -26.9     53.7      274 
## 
## Random effects:
## 
## Conditional model:
##  Groups  Name        Variance  Std.Dev. 
##  Site_ID (Intercept) 1.327e-09 3.643e-05
## Number of obs: 278, groups:  Site_ID, 16
## 
## Conditional model:
##                 Estimate Std. Error z value Pr(>|z|)    
## (Intercept)      -3.6689     0.9707  -3.780 0.000157 ***
## weigh_closeness  -2.3878    10.8177  -0.221 0.825305    
## niche.overlap     0.3487     3.7945   0.092 0.926773    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##  Family: poisson  ( log )
## Formula:          
## mean.seedweight ~ norm_degree + tot.visits.sc + weigh_closeness +  
##     niche.overlap + (1 | Site_ID)
## Data: s2
## 
##      AIC      BIC   logLik deviance df.resid 
##     65.4     87.2    -26.7     53.4      272 
## 
## Random effects:
## 
## Conditional model:
##  Groups  Name        Variance  Std.Dev. 
##  Site_ID (Intercept) 9.829e-10 3.135e-05
## Number of obs: 278, groups:  Site_ID, 16
## 
## Conditional model:
##                 Estimate Std. Error z value Pr(>|z|)    
## (Intercept)      -3.8657     1.1590  -3.335 0.000852 ***
## norm_degree       1.5590     3.6584   0.426 0.669994    
## tot.visits.sc     0.6076     2.5717   0.236 0.813219    
## weigh_closeness  -4.3360    12.0297  -0.360 0.718518    
## niche.overlap     0.1473     4.6368   0.032 0.974665    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##          dAIC df
## m2.seedw 0.0  4 
## m1.seedw 0.1  4 
## m3.seedw 3.7  6

##  Family: poisson  ( log )
## Formula:          
## mean.fruitweight ~ norm_degree + tot.visits.sc + (1 | Site_ID)
## Data: s2
## 
##      AIC      BIC   logLik deviance df.resid 
##    143.5    158.1    -67.8    135.5      277 
## 
## Random effects:
## 
## Conditional model:
##  Groups  Name        Variance  Std.Dev. 
##  Site_ID (Intercept) 9.866e-10 3.141e-05
## Number of obs: 281, groups:  Site_ID, 16
## 
## Conditional model:
##               Estimate Std. Error z value Pr(>|z|)    
## (Intercept)    -2.6777     0.4184  -6.400 1.55e-10 ***
## norm_degree     0.1324     1.9797   0.067    0.947    
## tot.visits.sc   0.5485     1.1716   0.468    0.640    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##  Family: poisson  ( log )
## Formula:          
## mean.fruitweight ~ weigh_closeness + niche.overlap + (1 | Site_ID)
## Data: s2
## 
##      AIC      BIC   logLik deviance df.resid 
##    142.1    156.6    -67.0    134.1      274 
## 
## Random effects:
## 
## Conditional model:
##  Groups  Name        Variance  Std.Dev.
##  Site_ID (Intercept) 1.005e-09 3.17e-05
## Number of obs: 278, groups:  Site_ID, 16
## 
## Conditional model:
##                 Estimate Std. Error z value Pr(>|z|)    
## (Intercept)      -2.4495     0.5347  -4.581 4.62e-06 ***
## weigh_closeness  -3.5110     6.0510  -0.580    0.562    
## niche.overlap     0.7745     2.0457   0.379    0.705    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##  Family: poisson  ( log )
## Formula:          
## mean.fruitweight ~ norm_degree + tot.visits.sc + weigh_closeness +  
##     niche.overlap + (1 | Site_ID)
## Data: s2
## 
##      AIC      BIC   logLik deviance df.resid 
##    145.8    167.5    -66.9    133.8      272 
## 
## Random effects:
## 
## Conditional model:
##  Groups  Name        Variance  Std.Dev. 
##  Site_ID (Intercept) 1.624e-09 4.029e-05
## Number of obs: 278, groups:  Site_ID, 16
## 
## Conditional model:
##                 Estimate Std. Error z value Pr(>|z|)    
## (Intercept)      -2.5218     0.6199  -4.068 4.74e-05 ***
## norm_degree       0.6631     2.1299   0.311    0.756    
## tot.visits.sc     0.4922     1.4143   0.348    0.728    
## weigh_closeness  -4.4036     6.6338  -0.664    0.507    
## niche.overlap     0.4962     2.4762   0.200    0.841    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##           dAIC df
## m2.fruitw 0.0  4 
## m1.fruitw 1.4  4 
## m3.fruitw 3.7  6

Network-level

FRUITSET

We find that at the level of the whole communty, differences in average fruit set are driven particularly by the diversity of pollinator species. In this sense, we find that communities with greater pollinator diversity have lower values of fruit set but consistent higher values of other function aspects measured: seed number per fruit, fruit and seed weight. In contrast, we found that metrics accounting for the structure of interactions between plants and pollinators had little or no effect for functioning.

##  Family: gaussian  ( identity )
## Formula:          mean.fruitset ~ tot.visits + species.poll
## Data: f.agg
## 
##      AIC      BIC   logLik deviance df.resid 
##      -22      -19       15      -30       12 
## 
## 
## Dispersion estimate for gaussian family (sigma^2): 0.00895 
## 
## Conditional model:
##                Estimate Std. Error z value Pr(>|z|)    
## (Intercept)   1.045e+00  9.874e-02  10.586  < 2e-16 ***
## tot.visits    9.685e-06  1.886e-05   0.513 0.607662    
## species.poll -6.267e-03  1.899e-03  -3.300 0.000966 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##  Family: gaussian  ( identity )
## Formula:          mean.fruitset ~ nodf.song + functional.comp.poll
## Data: f.agg
## 
##      AIC      BIC   logLik deviance df.resid 
##    -18.6    -15.5     13.3    -26.6       12 
## 
## 
## Dispersion estimate for gaussian family (sigma^2): 0.0111 
## 
## Conditional model:
##                      Estimate Std. Error z value Pr(>|z|)    
## (Intercept)          0.953432   0.124294   7.671 1.71e-14 ***
## nodf.song            0.900098   0.350255   2.570   0.0102 *  
## functional.comp.poll 0.001620   0.002403   0.674   0.5002    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##  Family: gaussian  ( identity )
## Formula:          
## mean.fruitset ~ tot.visits + species.poll + nodf.song + functional.comp.poll
## Data: f.agg
## 
##      AIC      BIC   logLik deviance df.resid 
##    -21.8    -17.2     16.9    -33.8       10 
## 
## 
## Dispersion estimate for gaussian family (sigma^2): 0.00707 
## 
## Conditional model:
##                        Estimate Std. Error z value Pr(>|z|)    
## (Intercept)           9.862e-01  1.014e-01   9.729  < 2e-16 ***
## tot.visits           -7.575e-06  1.876e-05  -0.404  0.68639    
## species.poll         -8.188e-03  2.717e-03  -3.014  0.00258 ** 
## nodf.song             2.607e-01  3.655e-01   0.713  0.47567    
## functional.comp.poll  4.735e-03  2.339e-03   2.025  0.04292 *  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##       dAIC df
## m1.fr 0.0  4 
## m3.fr 0.2  6 
## m2.fr 3.5  4

##  Family: nbinom2  ( log )
## Formula:          mean.seedset ~ tot.visits + species.poll
## Data: s.comm
## 
##      AIC      BIC   logLik deviance df.resid 
##    148.4    151.5    -70.2    140.4       12 
## 
## 
## Overdispersion parameter for nbinom2 family (): 1.69 
## 
## Conditional model:
##                Estimate Std. Error z value Pr(>|z|)    
## (Intercept)  -0.4697399  1.1219864  -0.419   0.6755    
## tot.visits    0.0008381  0.0002111   3.970  7.2e-05 ***
## species.poll  0.0701465  0.0214970   3.263   0.0011 ** 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##  Family: nbinom2  ( log )
## Formula:          mean.seedset ~ nodf.song + functional.comp.poll
## Data: s.comm
## 
##      AIC      BIC   logLik deviance df.resid 
##    154.1    157.2    -73.0    146.1       12 
## 
## 
## Overdispersion parameter for nbinom2 family (): 1.24 
## 
## Conditional model:
##                      Estimate Std. Error z value Pr(>|z|)  
## (Intercept)           1.01550    1.44933   0.701   0.4835  
## nodf.song            -0.50270    3.83045  -0.131   0.8956  
## functional.comp.poll  0.04554    0.02174   2.095   0.0362 *
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##  Family: nbinom2  ( log )
## Formula:          
## mean.seedset ~ tot.visits + species.poll + nodf.song + functional.comp.poll
## Data: s.comm
## 
##      AIC      BIC   logLik deviance df.resid 
##    151.7    156.4    -69.9    139.7       10 
## 
## 
## Overdispersion parameter for nbinom2 family (): 1.75 
## 
## Conditional model:
##                        Estimate Std. Error z value Pr(>|z|)   
## (Intercept)          -0.8786169  1.3199938  -0.666  0.50565   
## tot.visits            0.0007591  0.0002498   3.039  0.00237 **
## species.poll          0.0557510  0.0283785   1.964  0.04947 * 
## nodf.song            -0.7208715  4.1773613  -0.173  0.86299   
## functional.comp.poll  0.0173348  0.0242116   0.716  0.47401   
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##             dAIC df
## m1.seedcomm 0.0  4 
## m3.seedcomm 3.3  6 
## m2.seedcomm 5.7  4

##  Family: gaussian  ( identity )
## Formula:          mean.fruitweight ~ tot.visits + species.poll
## Data: s.comm
## 
##      AIC      BIC   logLik deviance df.resid 
##    -61.6    -58.5     34.8    -69.6       12 
## 
## 
## Dispersion estimate for gaussian family (sigma^2): 0.000755 
## 
## Conditional model:
##                Estimate Std. Error z value Pr(>|z|)   
## (Intercept)  -1.549e-02  2.920e-02  -0.531  0.59573   
## tot.visits    1.476e-05  9.662e-06   1.528  0.12655   
## species.poll  1.583e-03  5.509e-04   2.873  0.00407 **
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##  Family: gaussian  ( identity )
## Formula:          mean.fruitweight ~ nodf.song + functional.comp.poll
## Data: s.comm
## 
##      AIC      BIC   logLik deviance df.resid 
##    -58.7    -55.6     33.4    -66.7       12 
## 
## 
## Dispersion estimate for gaussian family (sigma^2): 0.000905 
## 
## Conditional model:
##                        Estimate Std. Error z value Pr(>|z|)
## (Intercept)          -0.0068309  0.0354505  -0.193    0.847
## nodf.song            -0.0731545  0.0998975  -0.732    0.464
## functional.comp.poll  0.0010219  0.0006853   1.491    0.136
##  Family: gaussian  ( identity )
## Formula:          
## mean.fruitweight ~ tot.visits + species.poll + nodf.song + functional.comp.poll
## Data: s.comm
## 
##      AIC      BIC   logLik deviance df.resid 
##    -57.8    -53.1     34.9    -69.8       10 
## 
## 
## Dispersion estimate for gaussian family (sigma^2): 0.000748 
## 
## Conditional model:
##                        Estimate Std. Error z value Pr(>|z|)
## (Intercept)          -2.113e-02  3.344e-02  -0.632    0.527
## tot.visits            1.360e-05  1.067e-05   1.274    0.203
## species.poll          1.361e-03  8.784e-04   1.550    0.121
## nodf.song            -1.189e-02  1.193e-01  -0.100    0.921
## functional.comp.poll  2.523e-04  7.522e-04   0.336    0.737
##           dAIC df
## m1.fwcomm 0.0  4 
## m2.fwcomm 2.9  4 
## m3.fwcomm 3.9  6

##  Family: gaussian  ( identity )
## Formula:          mean.seedweight ~ tot.visits + species.poll
## Data: s.comm
## 
##      AIC      BIC   logLik deviance df.resid 
##    -99.1    -96.0     53.5   -107.1       12 
## 
## 
## Dispersion estimate for gaussian family (sigma^2): 7.25e-05 
## 
## Conditional model:
##                Estimate Std. Error z value Pr(>|z|)    
## (Intercept)  -1.203e-02  9.051e-03  -1.329 0.183742    
## tot.visits    5.806e-06  2.995e-06   1.939 0.052537 .  
## species.poll  6.028e-04  1.708e-04   3.530 0.000415 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##  Family: gaussian  ( identity )
## Formula:          mean.seedweight ~ nodf.song + functional.comp.poll
## Data: s.comm
## 
##      AIC      BIC   logLik deviance df.resid 
##    -91.9    -88.8     49.9    -99.9       12 
## 
## 
## Dispersion estimate for gaussian family (sigma^2): 0.000114 
## 
## Conditional model:
##                        Estimate Std. Error z value Pr(>|z|)
## (Intercept)          -0.0007202  0.0125778  -0.057    0.954
## nodf.song            -0.0289963  0.0354437  -0.818    0.413
## functional.comp.poll  0.0002320  0.0002431   0.954    0.340
##  Family: gaussian  ( identity )
## Formula:          
## mean.seedweight ~ tot.visits + species.poll + nodf.song + functional.comp.poll
## Data: s.comm
## 
##      AIC      BIC   logLik deviance df.resid 
##    -95.8    -91.2     53.9   -107.8       10 
## 
## 
## Dispersion estimate for gaussian family (sigma^2): 6.93e-05 
## 
## Conditional model:
##                        Estimate Std. Error z value Pr(>|z|)   
## (Intercept)          -7.986e-03  1.018e-02  -0.785  0.43266   
## tot.visits            6.608e-06  3.249e-06   2.034  0.04195 * 
## species.poll          7.622e-04  2.674e-04   2.851  0.00436 **
## nodf.song             9.280e-03  3.632e-02   0.256  0.79834   
## functional.comp.poll -1.765e-04  2.289e-04  -0.771  0.44065   
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##           dAIC df
## m1.swcomm 0.0  4 
## m3.swcomm 3.3  6 
## m2.swcomm 7.2  4

EQUITY IN FRUIT PRODUCTION

LOOK AT EVENNESSS IN FUNCTION. within a site most plant species have very low reproductive success and thus we get incredibly high evenness values. therefore we use thresholds, for each sps and site we calculate the # of individuals that are above 50% in fruit production.

AS before we find that the best model is that including the diversity of pollinator species, which has a significant positive effect on equity in fruitset.

##  Family: gaussian  ( identity )
## Formula:          no_rows ~ tot.visits + species.poll
## Data: f6
## 
##      AIC      BIC   logLik deviance df.resid 
##     56.9     60.0    -24.5     48.9       12 
## 
## 
## Dispersion estimate for gaussian family (sigma^2): 1.25 
## 
## Conditional model:
##               Estimate Std. Error z value Pr(>|z|)    
## (Intercept)  0.7555909  1.1656620   0.648 0.516851    
## tot.visits   0.0007352  0.0002227   3.301 0.000962 ***
## species.poll 0.0500192  0.0224169   2.231 0.025660 *  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##  Family: gaussian  ( identity )
## Formula:          no_rows ~ nodf.song + functional.comp.poll
## Data: f6
## 
##      AIC      BIC   logLik deviance df.resid 
##     59.0     62.1    -25.5     51.0       12 
## 
## 
## Dispersion estimate for gaussian family (sigma^2): 1.42 
## 
## Conditional model:
##                      Estimate Std. Error z value Pr(>|z|)   
## (Intercept)           1.02950    1.40211   0.734  0.46280   
## nodf.song             5.24030    3.95107   1.326  0.18474   
## functional.comp.poll  0.08446    0.02710   3.116  0.00183 **
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##  Family: gaussian  ( identity )
## Formula:          
## no_rows ~ tot.visits + species.poll + nodf.song + functional.comp.poll
## Data: f6
## 
##      AIC      BIC   logLik deviance df.resid 
##     57.6     62.2    -22.8     45.6       10 
## 
## 
## Dispersion estimate for gaussian family (sigma^2): 1.01 
## 
## Conditional model:
##                       Estimate Std. Error z value Pr(>|z|)  
## (Intercept)          0.3848398  1.2134304   0.317   0.7511  
## tot.visits           0.0005429  0.0002246   2.417   0.0156 *
## species.poll         0.0414381  0.0325201   1.274   0.2026  
## nodf.song            4.9563131  4.3750998   1.133   0.2573  
## functional.comp.poll 0.0479159  0.0279978   1.711   0.0870 .
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##       dAIC df
## m1.eq 0.0  4 
## m3.eq 0.7  6 
## m2.eq 2.0  4

Now look at motifs and indirect interactions. To this end, we compare the roles of plant species at each site with low (<0.5) and high (>0.5) fruit set. We calculate for each plant species at each site the sum-normalised roles of all plant species in motifs up to 5 nodes and plot them using NMDS. ALthough we find a large overlap between the roles of both types ofplant species, we do find that plant species with high fruit set tend to be located in the higher values of NMDS1 which is related to network positions in which a generalist plant species interacts with highly specialized pollinator partners (e.g, positions 16 and 46, see Fig. 1 in bmotif paper, Simmons et al preprint)

## Run 0 stress 0.0579615 
## Run 1 stress 0.1186405 
## Run 2 stress 0.05796525 
## ... Procrustes: rmse 0.00142633  max resid 0.008008636 
## ... Similar to previous best
## Run 3 stress 0.05795964 
## ... New best solution
## ... Procrustes: rmse 0.001380013  max resid 0.007994088 
## ... Similar to previous best
## Run 4 stress 0.07238514 
## Run 5 stress 0.05796058 
## ... Procrustes: rmse 0.001283616  max resid 0.007866695 
## ... Similar to previous best
## Run 6 stress 0.05795968 
## ... Procrustes: rmse 1.458356e-05  max resid 5.477475e-05 
## ... Similar to previous best
## Run 7 stress 0.05796605 
## ... Procrustes: rmse 0.0003604755  max resid 0.0020523 
## ... Similar to previous best
## Run 8 stress 0.1307137 
## Run 9 stress 0.05795965 
## ... Procrustes: rmse 5.224754e-05  max resid 0.0003140176 
## ... Similar to previous best
## Run 10 stress 0.05796522 
## ... Procrustes: rmse 0.0003563191  max resid 0.002041866 
## ... Similar to previous best
## Run 11 stress 0.0708063 
## Run 12 stress 0.07080058 
## Run 13 stress 0.05796515 
## ... Procrustes: rmse 0.0003460244  max resid 0.001990277 
## ... Similar to previous best
## Run 14 stress 0.08906242 
## Run 15 stress 0.05796089 
## ... Procrustes: rmse 0.001227501  max resid 0.007630813 
## ... Similar to previous best
## Run 16 stress 0.07707273 
## Run 17 stress 0.1275027 
## Run 18 stress 0.05796563 
## ... Procrustes: rmse 0.0003840991  max resid 0.002093695 
## ... Similar to previous best
## Run 19 stress 0.1186401 
## Run 20 stress 0.07080604 
## *** Solution reached

## 
## Call:
## adonis(formula = mp2[, 3:48] ~ mean.fruitset, data = mp2, method = "bray") 
## 
## Permutation: free
## Number of permutations: 999
## 
## Terms added sequentially (first to last)
## 
##               Df SumsOfSqs MeanSqs F.Model      R2 Pr(>F)
## mean.fruitset  1    0.1558 0.15575  1.0844 0.01747  0.343
## Residuals     61    8.7615 0.14363         0.98253       
## Total         62    8.9172                 1.00000

Albrecht, Matthias, Bernhard Schmid, Yann Hautier, and Christine B Müller. 2012. “Diverse pollinator communities enhance plant reproductive success.” Proceedings of the Royal Society B: Biological Sciences 279 (1748): 4845 LP–4852. http://rspb.royalsocietypublishing.org/content/279/1748/4845.abstract.

Fründ, Jochen, Carsten F Dormann, Andrea Holzschuh, and Teja Tscharntke. 2013. “Bee diversity effects on pollination depend on functional complementarity and niche shifts.” Ecology 94 (9). Ecological Society of America: 2042–54. doi:10.1890/12-1620.1.

Potts, Simon G., Jacobus Christiaan Biesmeijer, Claire Kremen, Peter Neumann, Oliver Schweiger, and William E Kunin. 2010. “Global pollinator declines: trends, impacts and drivers.” Trends Ecol Evol 25: 345–53.